Search for: All records

Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (10^8.0M_⊙≤M_⋆ ≲ 10^12.0M_⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 10^10.0–10^11.0M_⊙, while galaxies with stellar mass < 10¹⁰M_⊙and >10¹¹M_⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. 

ABSTRACT The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (MZ*R) in the hydrodynamic simulations Illustris, TNG, and EAGLE (Evolution and Assembly of GaLaxies and their Environment) to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the MZ*R across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g. SFR, gas fraction, and atomic gas mass). Since we find this correlation exists in all three simulations, each employing a subgrid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ‘burstiness’ of feedback within galaxies. 

ABSTRACT We present radial gas-phase metallicity profiles, gradients, and break radii at redshift z = 0–3 from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, Rbreak, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of ∼2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that Rbreak occurs where the gas mixing time-scale is ∼10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic discs can be thought of in two-parts: a star-forming inner disc with a steep gradient and a mixing-dominated outer disc with a flat gradient, with the break radius marking the region of transition between them. 

Abstract Current methods of identifying the ionizing source of nebular emission in galaxies are well defined for the era of single-fiber spectroscopy, but still struggle to differentiate the complex and overlapping ionization sources in some galaxies. With the advent of integral field spectroscopy, the limits of these previous classification schemes are more apparent. We propose a new method for distinguishing the ionizing source in resolved galaxy spectra by use of a multidimensional diagnostic diagram that compares emission-line ratios with velocity dispersion on a spaxel-by-spaxel basis within a galaxy. This new method is tested using the Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph Galaxy Survey (SAMI) Data Release 3 (DR3), which contains 3068 galaxies atz< 0.12. Our results are released as ionization maps available alongside the SAMI DR3 public data. Our method accounts for a more diverse range of ionization sources than the standard suite of emission-line diagnostics; we find 1433 galaxies with a significant contribution from non-star-forming ionization using our improved method as compared to 316 galaxies identified using only emission-line ratio diagnostics. Within these galaxies, we further identify 886 galaxies hosting unique signatures inconsistent with standard ionization by Hiiregions, active galactic nuclei, or shocks. These galaxies span a wide range of masses and morphological types and comprise a sizable portion of the galaxies used in our sample. With our revised method, we show that emission-line diagnostics alone do not adequately differentiate the multiple ways to ionize gas within a galaxy. 

Abstract We present spectroscopic confirmation of candidate strong gravitational lenses using the Keck Observatory and Very Large Telescope as part of our ASTRO 3D Galaxy Evolution with Lenses ( AGEL ) survey. We confirm that (1) search methods using convolutional neural networks (CNNs) with visual inspection successfully identify strong gravitational lenses and (2) the lenses are at higher redshifts relative to existing surveys due to the combination of deeper and higher-resolution imaging from DECam and spectroscopy spanning optical to near-infrared wavelengths. We measure 104 redshifts in 77 systems selected from a catalog in the DES and DECaLS imaging fields ( r ≤ 22 mag). Combining our results with published redshifts, we present redshifts for 68 lenses and establish that CNN-based searches are highly effective for use in future imaging surveys with a success rate of at least 88% (defined as 68/77). We report 53 strong lenses with spectroscopic redshifts for both the deflector and source ( z src > z defl ), and 15 lenses with a spectroscopic redshift for either the deflector ( z defl > 0.21) or source ( z src ≥ 1.34). For the 68 lenses, the deflectors and sources have average redshifts and standard deviations of 0.58 ± 0.14 and 1.92 ± 0.59 respectively, and corresponding redshift ranges of z defl = 0.21–0.89 and z src = 0.88–3.55. The AGEL systems include 41 deflectors at z defl ≥ 0.5 that are ideal for follow-up studies to track how mass density profiles evolve with redshift. Our goal with AGEL is to spectroscopically confirm ∼100 strong gravitational lenses that can be observed from both hemispheres throughout the year. The AGEL survey is a resource for refining automated all-sky searches and addressing a range of questions in astrophysics and cosmology. 

ABSTRACT We present the radial gas-phase, mass-weighted metallicity profiles and gradients of the TNG50 star-forming galaxy population measured at redshifts z = 0–3. We investigate the redshift evolution of gradients and examine relations between gradient (negative) steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity among these negative gradients. We determine that the gradients of all galaxies grow more negative with redshift at a roughly constant rate of approximately $$-0.02\ \mathrm{dex\, kpc^{-1}}/\Delta z$$. This rate does not vary significantly with galaxy mass. We observe a weak negative correlation between gradient (negative) steepness and galaxy stellar mass at z < 2. However, when we normalize gradients by a characteristic radius defined by the galactic star formation distribution, we find that these normalized gradients do not vary significantly with either stellar mass or redshift. We place our results in the context of previous simulations and show that TNG50 high-redshift gradients are more negative than those of models featuring burstier feedback, which may further highlight high-redshift gradients as important discriminators of galaxy formation models. We also find that z = 0 and z = 0.5 TNG50 gradients are consistent with the gradients observed in galaxies at these redshifts, although the preference for flat gradients observed in redshift z ≳ 1 galaxies is not present in TNG50. If future JWST (James Webb Space Telescope) and ELT (Extremely Large Telescope) observations validate these flat gradients, it may indicate a need for simulation models to implement more powerful radial gas mixing within the ISM (interstellar medium), possibly via turbulence and/or stronger winds. 

Abstract The 3D geometries of high-redshift galaxies remain poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in James Webb Space Telescope Cosmic Evolution Early Release Science observations with $\log M_{*}/M_{\odot }=9.0–10.5$atz= 0.5–8.0. We reproduce previous results from the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size, and covariances with samples as small as ∼50 galaxies. We find high 3D ellipticities for all mass–redshift bins, suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be ∼1 for $\log M_{*}/M_{\odot }=9.0–9.5$dwarfs atz> 1 (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a “banana” in the projected $b/a–\log a$diagram with an excess of low-b/a, large- $\log a$galaxies. The dwarf prolate fraction rises from ∼25% atz= 0.5–1.0 to ∼50%–80% atz= 3–8. Our results imply a second kind of disk settling from oval (triaxial) to more circular (axisymmetric) shapes with time. We simultaneously constrain the 3D size–mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices (n∼ 1), nonparametric morphological properties, and specific star formation rates. Both tend to be visually classified as disks or irregular, but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects, and theoretical implications. 

ABSTRACT We investigate the mean locally measured velocity dispersions of ionized gas (σgas) and stars (σ*) for 1090 galaxies with stellar masses $$\log \, (M_{\!\ast }/M_{\odot }) \ge 9.5$$ from the SAMI Galaxy Survey. For star-forming galaxies, σ* tends to be larger than σgas, suggesting that stars are in general dynamically hotter than the ionized gas (asymmetric drift). The difference between σgas and σ* (Δσ) correlates with various galaxy properties. We establish that the strongest correlation of Δσ is with beam smearing, which inflates σgas more than σ*, introducing a dependence of Δσ on both the effective radius relative to the point spread function and velocity gradients. The second strongest correlation is with the contribution of active galactic nuclei (AGN) (or evolved stars) to the ionized gas emission, implying that the gas velocity dispersion is strongly affected by the power source. In contrast, using the velocity dispersion measured from integrated spectra (σap) results in less correlation between the aperture-based Δσ (Δσap) and the power source. This suggests that the AGN (or old stars) dynamically heat the gas without causing significant deviations from dynamical equilibrium. Although the variation of Δσap is much smaller than that of Δσ, a correlation between Δσap and gas velocity gradient is still detected, implying that there is a small bias in dynamical masses derived from stellar and ionized gas velocity dispersions. 

ABSTRACT We measure the gas-phase metallicity gradients of 248 galaxies selected from Data Release 2 of the SAMI Galaxy Survey. We demonstrate that there are large systematic discrepancies between the metallicity gradients derived using common strong emission line metallicity diagnostics. We determine which pairs of diagnostics have Spearman’s rank coefficients greater than 0.6 and provide linear conversions to allow the accurate comparison of metallicity gradients derived using different strong emission line diagnostics. For galaxies within the mass range 8.5 < log (M/M⊙) < 11.0, we find discrepancies of up to 0.11 dex/Re between seven popular diagnostics in the metallicity gradient–mass relation. We find a suggestion of a break in the metallicity gradient–mass relation, where the slope shifts from negative to positive, occurs between 9.5 < log (M/M⊙) < 10.5 for the seven chosen diagnostics. Applying our conversions to the metallicity gradient–mass relation, we reduce the maximum dispersion from 0.11 dex/Re to 0.02 dex/Re. These conversions provide the most accurate method of converting metallicity gradients when key emission lines are unavailable. We find that diagnostics that share common sets of emission line ratios agree best, and that diagnostics calibrated through the electron temperature provide more consistent results compared to those calibrated through photoionization models. 

Abstract We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission-line maps at ∼25 pc resolution from the Hubble Space Telescope, Keck/NIRC2 with Adaptive Optics, and the Atacama Large Millimeter/submillimeter Array (ALMA). NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and filaments seen in shock tracers from warm molecular gas (H₂2.12μm) to optical ionized gas ([Oiii], [Nii], [Sii], and [Oi]) and hot plasma (FeXXV). In the most distinct bubble, we see a clear shock front traced by high [Oiii]/Hβand [Oiii]/[Oi]. Cool molecular gas (CO(2−1)) is only present near the base of the bubble, toward the nuclei launching the outflow. We interpret the lack of molecular gas outside the bubble to mean that the shock front is not responsible for dissociating molecular gas, and conclude that the molecular clouds are partly shielded and either entrained briefly in the outflow, or left undisturbed while the hot wind flows around them. Elsewhere in the galaxy, shock-excited H₂extends at least ∼4 kpc from the nuclei, tracing molecular gas even warmer than that between the nuclei, where the two galaxies’ interstellar media are colliding. A ridgeline of high [Oiii]/Hβemission along the eastern arm aligns with the southern nucleus’ stellar disk minor axis; optical integral field spectroscopy from WiFeS suggests this highly ionized gas is centered at systemic velocity and likely photoionized by direct line of sight to the southern active galactic nucleus.